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Abstract 



Superkamiokande atmospheric neutrino data (contained and upward going 
through [i events) for 990 days are analyzed in the framework of four neutri- 
nos (three active and one sterile neutrinos) without imposing constraints of 
£T) ■ Big Bang Nucleosynthesis. It is shown that the wide range of the oscillation 

parameters is allowed at 90% confidence level (0.1 < \U s i\ 2 + \U S 2\ 2 < 1), 
\ where the best fit point has some contribution of <-> v s and contribution 

of AmL SND . The case of pure «-> v s oscillation is excluded at 99.7%CL 
' (3. Oct) which is consistent with the recent analysis by the Superkamiokande 

group. Combining this result with the analysis by Giunti, Gonzarez-Garcia 



Ph. and Peha-Garay, it is found that the Large Mixing Angle and Vacuum Oscil- 

r~| ! lation solutions of the solar neutrino problem are also allowed. 
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I. INTRODUCTION 



There have been a few experiments which suggest neutrino oscillations: the solar neutrino 
deficit the atmospheric neutrino anomaly |7|-[T6|] and the LSND data [17]. If we assume 
that all these three are caused by neutrino oscillations then we need at least four species 
of neutrinos and schemes with sterile neutrinos have been studied by many people |T8|-|25 



By generalizing the discussion on Big Bang Nucleosynthesis (BBN) from the two neutrino 
scheme f26| to four neutrino case, it has been shown |2(J [jthat the neutrino mixing angles are 
strongly constrained not only by the reactor data |T^] but also by BBN if one demands that 



the number N u of effective neutrinos be less than four. In this case the 4x4 MNS matrix 27 



splits approximately into two 2x2 block diagonal matrices, and the solar neutrino deficit 
is explained by u e <-> v s oscillations with the Small Mixing Angle (SMA) MSW solution 
|28| and the atmospheric neutrino anomaly is accounted for by <-> u T . On the other 
hand, some people have given conservative estimate for N v |29| and if their estimate is 
correct then we no longer have strong constraints on the mixing angles of sterile neutrinos. 
Recently Giunti, Gonzalez-Garcia and Pena-Garay |3(| have analyzed the solar neutrino data 
in the four neutrino scheme without BBN constraints. They have shown that the scheme is 
reduced to the two neutrino framework in which only one free parameter c s = \U s i\ 2 + \ U S 2\ 2 
appears in the analysis. Their conclusion is that the SMA MSW solution exists for the 
entire region of < c s < 1, while the Large Mixing Angle (LMA) and Vacuum Oscillation 
(VO) solutions survive only for < c s ^ 0.2 and < c s ^ 0.4, respectively. In this paper 
the Superkamiokande atmospheric neutrino data (contained and upward going through fi 
events) are analyzed in the same scheme as in ||30|| , i.e., in the four neutrino scheme with 
all the constraints of accelerators and reactors and without BBN constraints. Zenith angle 
dependence of atmospheric neutrinos has been analyzed by many theorists pT] , |23| -p5|] as well 
as by experimentalists [M||[l0Jl3]-O . 



II. FOUR NEUTRINO SCHEME 

Here we adopt the notation in [[2(J for the 4x4 MNS matrix: 











= U 














w 



(Vel U e2 U e3 U e ,\ 

Ufj,l U^2 U^,3 

U T 1 U T 2 U T s U T 4 

\u sl u s2 u s3 U S J 



= R 3 a(~ - e 34 )R 24 (e 24 )R 23 ( - ) e 2i5lA 3 J R 23 (^ 3 )e- 2i5lA3 e^ Al5 / 2 

x i? 1 4(^i4)e-^ 3Al5/2 e 2i52A8/ ^ J Ri3(^i3)e- 2i52A8/v/5 i?i 2 (0i2), (1) 



x The result in [2C] was refined later in [22| with more careful treatment 



where Cy = cos = sin^j and 



i2j fc (0) = exp (iTj-a 



(2) 



is a 4 x 4 orthogonal matrix with 



(T 3 



km 



8 jm.3 kl) , 



(3) 



and 2A 3 = diag(l, -1, 0, 0), 2^3A 8 = diag(l, 1, -2, 0), 2y/6X 15 = diag(l, 1, 1, -3) are diago- 
nal elements of the 577(4) generators. 

We can assume without loss of generality that 



m\ < m 2 < m 2 < m\. 



(4) 



Three mass scales Am 2 ~ £>(l(r 5 eV 2 ) or e>(l(T 10 eV 2 ), Am 2 tm ~ C(10~ 2 eV 2 
(9(leV 2 ) are necessary to explain the suppression of the 7 Be solar neutrinos [ 
angle dependence of the atmospheric neutrino data fllO|, and the LSND data [17 



a ™lsnd ~ 
the zenith 



so we 
0T) are 



assume that three independent mass squared differences which are obtained from 
AitLq, Am 2 tm , Am 2 SND . It has been known [^,^l|] that schemes with three degenerate 
masses and one distinct massive state do not work to account for all the three neutrino 
anomalies, but schemes with two degenerate massive states (m 2 ~ m\ <C m 2 ~ m 2 , where 
(a) (Am 2 1; Am 2 3 ) = (Am|, Am 2 tm ) or (b) (Am 2 1; Am 2 3 ) = (Am 2 tm , Am|)) do. As far as 
the analyses of atmospheric neutrinos and solar neutrinos are concerned, the two cases (a) 
and (b) can be treated in the same manner, so we assume Am^ = Affig, Am 2 3 = Am 2 tm 
in the following. 

For the range of the Am 2 suggested by the LSND data, which is given by 0.2 eV 2 ^$ 
Am 2 SND < 2 eV 2 when combined with the data of Bugey and E776 |]3"5fl , the constraint 
by the Bugey data is very stringent and 



\U e3 \ 2 + \U e4 \ 2 <1U~ 2 , 



(5) 



has to be satisfied p9|-pT|, therefore we put U e % = U e ^ = for simplicity in the following 
discussions. Also in the analysis of atmospheric neutrinos, \Am 2 ? ) L / AE\ <d 1 is satisfied for 
typical values of the neutrino path length L and the neutrino energy E for atmospheric 
neutrinos, so we assume Amg = for simplicity throughout this paper. 

In the following analysis we will consider the situation where non-negligible contribu- 
tions from the largest mass squared difference Am 2 SND appear in the oscillation probability 
P{v„. — > V,,)- To avoid contradiction with the negative result of the CDHSW disappearing 

0.3eV 2 clS db reference value. 



'fj, 

experiment on v. 
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we will take Am 2 SND = 



Having assumed U e3 = U ei = and Am^ = 0, we have only mixings among u T , v s in 
the analysis of atmospheric neutrinos, and the Schrodinger equation we have to consider is 



d 

dx 



v T (x) 



f/diag (-AE 32 , 0, AE 43 ) U- 1 + diag (0, 0, A{x)) 



u T (x) 



(6) 



where 



3 



u 



Ut2 U t s Uri 

\u s2 u s3 u s4 



-C 2 4S 2 3e 



AS 



C23C24 



C23S34 + C 3 4S 2 3' S 24e _C 23C34S24 ~ S23S34C 



C23C34 + S23S24S34e il51 -C23S24S34 + C34S23C 



-iSi 
-iSi 



S24 
C24C34 
C24S34 



(7) 



AEij = Am-j-/2£ and A(x) = GpN n (x)/\^2 stands for the effect due to the neutral 



p6| after adding the unit 
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current interactions between z/ M , v T and matter in the Earth 
matrix diag (A(x), A(x), A(x)) to the right hand side of (H). Since v e does not oscillate 
with any other neutrinos, the only oscillation probability which is required in the analy- 



sis of atmospheric neutrinos is Piy^ 



In the present case, we have mass hierarchy 



|AE 32 | > |A£ 43 |, \A{x)\ and to the leading order in | A£ 43 |/| A£ 32 | and |A(x)|/|A£ 32 | we 
can obtain the analytical expression for the oscillation probability P{p i1 — > Vp) in adiabatic 
approximation, i.e., assuming that the derivative \dA(x)/dx\ is not large compared to AEjk- 
As we will see in the appendix, with such approximation we can show that the oscillation 



probability P(fu 



is invariant under # 34 — » — ^34 and 5\ — > n — 5\ to the leading order in 



|A_E 32 | and |A(x)|/|A_E 32 |. So we will study our scheme for the range < 8 2 4 < 7r/2, 
< 7r/2, — 7r/2 < # 34 < 7r/2, < 5\ < ix / 2. Note that all the three mixing angles 



|A£ 43 |/ 

< #23 

would lie in the first quadrant |35| if there were no matter effects. 



III. ANALYSIS OF THE ATMOSPHERIC NEUTRINO DATA 

We calculate the disappearance probability Piy^ — > u^) by solving numerically, and 
evaluate the number of events: 

N(n) = (l-r M )iV (/i)+r^iVo( e ), 
N(e) = (l-r e )N (e)+r e N (ji), 

N (ii) = n T J dE J dq J dcosQ J dcosO J dip e At (g) 

x d 3 F,(E,9) d 2 a,(E,q) dco^ ^ v ■ E 9) 

dE d cos 9 dip dq d cos ip d cos 9 ' ' ' 

N (e) = rix J dE J dq J dcosO J dcos9 J dp) e e (q) 

d 3 F e (E,9) d 2 a e (E,q) d cos ip 
dE d cos 9 d(p dq d cos ip d cos 9 ' 

where d 3 F a /dEd cos 9 'dp (a = /1, e) is the flux of atmospheric neutrinos u a (a = fi,e) with 
energy E from the zenith angle 9 nx is the effective number of target nucleons, e a (q) 
is the detection efficiency function for charged leptons a (a = /i, e) [1381 , da a jdqd cos ip 
(a = /i, e) is the differential cross section of the interaction v a N — > al (a = e or //; for sub- 
GeV events quasi-elastic scatterings z/ a iV — > aN' are dominant and the cross-section given 
in [3J| is used, while for multi-GeV events the inclusive cross-section for u a N — > aX given 
in HO] is used), and is the zenith angle of the direction from which the charged lepton a 
comes (See Fig. 1 of |P3j|), r M and r e stand for ratios of contamination due to misidentification 
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of /i-like and e-like events and we have put 7^=0.03, r e =0.06 for sub-GeV events, r A1 =0.007, 
r e =0.12 for fully contained multi-GeV events [11], 7^=0 for partially contained multi-GeV 
and upward going through \i events. 
X 2 is defined as 



2 _ 2 I 2 I 2 

X — Xsub-GeV + Amulti-GeV Xth_ 



rough 



(9) 



where 



Xsub-GeV 



Xmulti-GeV 



Xthrough 




1 - f )N](e) 



n°{ej 



+ 



a 



;i + f)iV|( M )-n5( M ) 



%0^f (e) - nf(t 



rij{e) 



a 



th 



10 



th 



+ E 



[a th iVj(/i) - n){y) 



3=1 



+ 



n§(//) 

a(l + ^)A^ m (/i) - nf (//) 



(10) 



are x 2 for sub-GeV, multi-GeV, and upward going through // events, respectively, the sum- 
mation on j runs over the ten bins for each % 2 , Nj(a) and (a=s, m, th) stand for 
theoretical predictions and data for the numbers of sub-GeV, multi-GeV, and upward going 
through n events, and it is understood that \ 2 is minimized with respect to all the nor- 
malization factors a, /3 S , f3 m , a t h- We have put a s =0.08, a m =0.12, cr ath =0.22 and we have 
assumed that the overall flux normalization a in the contained events is a free parameter as 
in ||10|| , and we have omitted the uncertainties of E u spectral index, relative normalization 
between PC and FC and up-down correlation for simplicity 

We have evaluated x 2 for 8 24 = (25 + 5j)° (j = 0, ■ • • , 7), 6 34 = 15j° (j = -6, • ■ ■ , 6), 
# 23 = \0j° (j = o, • • • , 4), Si = 0°, 45°, 90°, Am 2 3 = l(r 4+J '/ 10 eV 2 (j = 5, • • • , 20) and it is 
found that \ 2 has the minimum value 



v 2 . 

A mm 



43.1 



(Xsub-GeV ~~ 19.0, Xmulti-GeV ~~ 13.2, ^through 



11.6) 



(11) 



for 



Am 2 3 



10 -2.9 e y2 = L3 x 1Q -3 eV 2 



'24, #34, #23 ) = (35°, 15°, 20°) 



(12) 



for 45 degrees of freedom. For pure 



?34 = #23 = 0), the best fit is obtained 



Xminl^M 



48.3 (Xsub-GeV ~~ 19.8, Xmulti-GcV ~~ 17.0, x\a 



rough 



10.6) 



(13) 



for 



Am 2 3 



2.0 x 10 _d eV z 



?24, #34, #23; 



(40°,0°,0°). 



(14) 



The allowed regions at lcrCL, 90%CL, 99%CL are obtained by x 2 < X m in + ^X 2 , where 
Ax 2 =5.9, 9.2, 15.1 for five degrees of freedom, respectively, and they are depicted in Fig. 1 
for various values of 6> 24 and for 8 X = (1(a)), 8 X = 7r/4 (1(b)) and 8 X = 7r/2 (1(c)). At 99% 
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confidence level we find 27° ^ #24 ^ 58°. The difference of the best fit point among the all 
parameter space and the best fit for pure z/ M <-> v T (Am^ = 10~ 27 / 10 eV 2 , 6*24 = 40°, 6*34 = 0°, 
$23 — 0°) is A% 2 =5.1 which corresponds to 0.84a, so this difference is not significant. Zenith 
angle dependence for contained events and upward going through muon events for three sets 
of the parameters including the case ( |i~2|) and ( |T4"D is given in Fig. 2(a) and 2(b). Zenith 
angle dependence for no oscillation case that we obtained is reasonably in good agreement 
with the Superkamiokande result |T2| and this puts confidence to the present analysis. In 
general, the reason that the best fit point is slightly away from pure <-> v T case is because 
a better fit to the multi-GeV contained events compensates a worse fit to the upward going 
through fi events, and in total the case of hybrid oscillations of <-> v T and <-» v s fits 
better to the data (See Figs. 2(a) and 2(b)). 



IV. DISCUSSIONS 

(1) Pure Up <-> v s oscillation is obtained only for #23 = 0, #34 = ±90°, and in this case x 2 
satisfies x 2 ~ Xmm >17.9 for any value of Am^, # 2 4 and 8\ and it implies that pure z/ M <-> v s 
oscillation is excluded at 99.7%CL (3.0cxCL). This is consistent with the recent claim by 
Superkamiokande group 

(2) For some Am 2 3 , 9 24 5 X (e.g., Arr4=1.3xl0- 3 eV 2 , # 24 = 50°, 9 23 = 30°, 5i = 0), 
#34 = 90° is allowed at 90%CL, and one might wonder how such a solution can give a 
reasonable fit to the data. For this set of the oscillation parameters, we find that the 
oscillation probabilities without matter effects are given by 

3 ,/5t 

PyVn — > v T \ without matter effects) = - cos 



' 57T \ 

P(i' n — without matter effects) sin 2 ( —J 



1 3 9 /Am 2 L s 



(15) 



where (3/8) cos 2 (57r/18) ~ 0.15, sin 2 (57r/9) ~ 0.97, and we have averaged over rapid oscil- 
lations which come from Am 2 2 . From (|i~5f) we observe that zenith angle dependence with a 
non-maximal coefficient comes solely from <-> v s oscillation and both <->■ and v s 



oscillations have constant contribution in the oscillation probabilities. As was shown in |23 
the zenith angle dependence of contained events is explained well also by <-> v s oscillation, 
but the fit of «-> v s to upward going \i events is poor [p],^2|. In the present case, we have 
xLb-GeV^S-T 7 ) Xmuiti-Gev = 15.5, Xthrou g h = 15.0 and we observe that this solution is allowed 
at 90%CL because a poor fit to the upward going through // events is compensated by a 
good fit to the multi-GeV contained events. 

(3) One of the remarkable features of our result is that there is a region where relatively 
small value of #24 (~ 30°) is allowed at 90%CL for #23 — 20°, i.e., there exists a solution in 
which all the mixing angles are relatively small \Qjk\ ~ tt/6, and such situation does not 
occur for #23 = 0. This phenomena is reminiscent of the work [f|l] on three flavor analysis of 
the atmospheric neutrino data, in which it was argued that relatively small mixing angles 
could account for the data because of matter effects if there were no CHOOZ constraint 



42]. However, in the present case it turns out that the matter effect is not so important, 



i.e., even if we put the density of the matter to zero the fit to the data is still reasonable. 
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The reason that the fit to the data is good for 624 ~ 30° is because the disappearance 
probability behaves as 1 — Piy^ «-» v^) = A + B sin 2 (Am| 3 L/4£') (A, B are constant) and 
even if the coefficient B is relatively small the presence of the constant term A compensates 
the goodness of the fit to a certain extent. This argument applies both to the multi-GeV 
/i-like contained events and the upward going through muon events. 

(4) To combine the present result with the analysis of the solar neutrinos in ]30|], it is 



necessary to obtain the value of c s = \U s i\ 2 + \U S 2\ 2 for each point. In our parametrization 
(H) we have 

U 8 i = —£12(023034 + S23S24S34^ Sl ) 

U s2 = C1 2 (C23C 3 4 + S23S24S 3 4e lSl ) (16) 

so that 

c s = \U sl \ 2 + \U s2 \ 2 

= |C23C 3 4 + S23S24S 3 4e i<5l | 2 . (17) 

The contours of c s =0.2, 0.4, 0.6, 0.8 are plotted together with the allowed region for various 
Amf 3 in Fig. 3(a) (Si = 0), 3(b) (fa = it/A) and 3(c) (6t = tt/2). For each value of 8 X 
some point in the allowed region satisfies c s ^ 0.2 for 40° ^ 824 ^ 50°. We find that the 
allowed region at 90% confidence level satisfies 0.15 ^ c s < 1 for 5\ = 0, 0.10 ^ c s < 1 for 
5\ = 7r/4, 0.05 ^ c s < 1 for 5i = tt/2, respectively. Hence combination of the present result 
with the analysis in |30| suggests that the LMA and VO solutions as well as SMA solution 



of the solar neutrino problem are possible for some region in the parameter space. Recently 
Superkamiokande group has announced Jl3 that the SMA and VO solutions of the solar 



neutrino problem are disfavored at 95% confidence level. If Nature is described by a four 
neutrino scenario, therefore, the present scheme with c s ^ 0.2 may be the right solution for 
all the oscillation data. 



V. CONCLUSIONS 

We have shown in the framework of four neutrino oscillations without assuming the 
BBN constraints that the Superkamiokande atmospheric neutrino data are explained by 
wide range of the oscillation parameters which implies hybrid oscillations of <-» v T and 
v s as well as hybrid oscillations with Am 2 tm and AmL SND . The case of pure <-» v s 
is excluded at 3.0<rCL in good agreement with the Superkamiokande analysis. It is found 
by combining the analysis on the solar neutrino data by Giunti, Gonzalez-Garcia and Pena- 
Garay that the LMA and VO solutions as well as SMA solution of the solar neutrino problem 
are allowed. This gives us another possibility in phenomenology of neutrino oscillations and 
such scenarios deserve further study. 
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APPENDIX 



With hierarchy |Ai? 32 | ^> |A_E 43 |, |A(a;)| we can obtain the analytical expression for the 
oscillation probability Piy^ — > v^) to the leading order in | AE i3 \/\ AE 32 \ and |^4(x) |/| A_E 32 | 
in adiabatic approximation. The purpose of this Appendix is to show that the oscillation 
probability Piy^ <-> u^) is invariant under # 34 — > —#34 and 5i — > n — 5\ to the leading order in 
I AE 43 /AE 32 \ and \A{x) / AE 32 \. To do that, it is convenient to use different parametrization: 

V = R 13 (~) VR 23 (<p 23 )R 13 (<p 13 )V- 1 R 12 (i Pl2 ), (18) 

where V = diag(e*' 5 ',l,l). To relate this parametrization with (|7|), we have to multiply 
diagonal phase matrices 

V = diag (e lLU \e iuj \ e ~ i( " 1+W2) ) [/diag (V 71 , e* 72 , e - l(7l+72) ) , (19) 

where 

= - arg £/ s2 + ^ arg £/ s3 - - arg [7^, ^2 = ^ arg [7 s2 + ^ arg £/ s3 + arg 

12 1 2 11 

7i = -3 argf/ s2 + -argf/ s3 + -arg^, 72 = - arg C/ s2 - - arg U s3 + - arg U^. (20) 

Comparing the both hand sides of (|T9|), we obtain 

1 0.3 C24S34 

tan<£> 12 = , tan (^13 



,^2 1 ' ^Jc 2 2i c 2 u + s 2 2i 

tan v9 23 = , arg (-V e3 ) = it - 5' = - arg - arg U s2 - arg C/ s3 . (21) 

Using this parametrization, we have 

7 _1 diag (0, 0,A)V+ (0, 0, AE 43 ) 
= RnM~ l VR l3 { Vl3 - V^-Miag (*-, 0, *+) Rm&is - <f^)V- 1 R 12 ( (fl2 ) 
= diag (c 2 l2 (t^C 2 + t+5 2 ), S 1 2 2 (t_C' 2 + t + S 2 ), t+C 2 + i_S 2 ) 

+ ^sin2v9i 2 (t_C' 2 + t+S' 2 )Ai + -C 12 (t_ - t+) sin2<p(A 4 cos5' - A 5 sin5') 

+ -5'i2(i- — £+) sin2</5(A6 cos 5' — Aysin^'), (22) 

where Xj(j = 1, • • • , 8) are the 3x3 Gell-Mann matrices with normalization tr(AjAfc) = 25jk, 
S12 = sin<y9 12 , Cxi = cos<^i 2 , 5 ee sin(<^ 13 — y?ff), C* = cos(y?i 3 — <£>*|) and ip^ is the effective 
mixing angle in matter given by 

tan2^EE a ^ sin 2^3 

^ 13 AE 43 cos 2^9 13 — A K ' 

and t± are the eigenvalues defined by 
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t ± = -[A + AE A3 ± V(AS 43 cos 2<p 13 - A) 2 + (A£ 43 sin 2<^ 13 ) 2 ) . 



(24) 



Putting 



A = 



2AE 



32 



sin 2Lpi 2 (t-C 2 + t + S 2 )\i + Cu(t- — £+) sin2<^(A4Cos5' — A5 sin 5') 



(25) 



we have 



-iA 



diag (-A£ 32 , 0, 0) + 7 _1 diag (0, 0, A) V + diag (0, 0, AE 43 ] 



JA 



diag (-AE32, 0, 0) + (C 2 12 (t-C 2 + t+£ 2 )/ 3 



1 r 



{t_C 2 + t+S 2 )^ - 2C 1 2 2 ) + i+C* 2 + t_S 2 J diag(0, 1, 1) 

^(t-^ + t+S 2 ) -t + C 2 -t^S 2 ] \ 8 + ~S 12 (t_ - t + ) sm2<p(\ 6 cos 5' - \ 7 sin 5'), (26) 



+ 2 
1 

+ 2 



where J3 = diag(l, 1, 1), Ag = diag(0, 1, —1). The last two lines in ( p6|) can be diagonalized 
by noting 



JFA 8 + <?(A 6 cos 5' — A 7 sin 5') 



TXs + Qe 



A fi e 2 



<5'A 8 



As 



(27) 



where 



1 



sUt~c 2 + t + S 2 )-t + C 2 - t_S< 



Q = ^S 12 {t- - t + ) sin2<^ 

and the effective mixing angle ip in matter is given by 

Q 



tan 2ip 



(28) 



(29) 



Putting everything together, we obtain 

Vdiag (-AS32, 0, AE43) + diag (0, 0, A) 
= Ve ¥ r ^ e -^e iK {diag (-AS32, 0, 0) + (C 2 2 (t_C- 2 + t + S 2 )/ 3 



1 

+ 2 



(i_(7 2 + i+5 2 )(5 2 2 - 2C\ 2 ) + t+C* 2 + t_5 2 
xe -iA e i^x 7e -±8'\s V -i 



diag(0, 1, 1) + VF 2 + Q 2 \ 



(30) 



to first order in \AE i3 / AE 32 \ and \A(x)/ AE 32 \. The effective MNS matrix in matter is 
given by Ve lS ' Xs ^ 2 e~' i ^ X7 to zeroth order, and the three eigenvalues in matter are —AE 32 , 
f+, f_, where f± = (t_C* 2 + t+5 2 )(S 2 2 - 2C 2 2 ) + £ + C* 2 + t_S 2 ± V-? 72 + G 2 , to first order in 
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IAE43/A-E32I and \A{x)/ AE 32 \, after subtracting the contribution from Cf 2 {t~C 2 + t + S 2 )I 3 . 
The oscillation probability P(f M — > z^) is given by 



P(j^n - = 1 - 4|^| 2 (1 - |^| 2 ) sin 2 

-4|e^'CVV; 2 + e^'S^H - e^S^, + e -**'C^| 2 sin 2 



(31) 

to the leading order in | A.E 43 |/| AE 32 \ and |A(a;)|/|A£/32|, where = cosip, = smtp and 
V^j are the matrix elements given by (|IBD : 

V^i = — S12C23 — Ci2S23Sise lS 

Vfj,2 = C12C23 — 5'i25'23S'i3e" 5 

Vn3 = S23C13 (32) 

with = sinyjjfc, (T^ = cos^jfc- We observe that (|3l| ) is invariant under #34 — > —6*34 and 
<5i — 7r — <5i since under these transformation the mixing angles in V behave as tp\2 — > tpu, 
V13 -> -Via (^13 -> -^13), ^23 -> V23, argU sj -> - arg C/ Sj (j=2,3), argC^ 4 -> tt - argC/^, 
5' — 7T - 5', -0 -> ^ - V; 3 (.7 = 1,2), V^3 -> \ei 5 'C^ 2 + e-^'S^l 2 -> 

le^'C^ + I - e^'^^2 + e^'C^I 2 - I - e^S^ + e~^' C^V^ 2 ■ 

Therefore, the oscillation probability is invariant under #34 — > —#34, 5i — > 7r — 5i to zeroth 
order in \AE i3 /AE 32 \ and \A(x)/AE 32 \. 
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Figures 



Fig.l (a), (b), (c) The allowed region in the (#34, #23) plane for various values of #24 = 
27°,30V--,58° and (a) S ± = (Am 2 3 = l(r 29 / 10 eV 2 ), (b) ^ = tt/4 (Am 2 3 = 
10- 27 /io e v 2 ), (c) 6 1 = tt/2 (Am 2 3 = 10- 29 / 10 eV 2 ), respectively. The solid, dashed 
and dotted lines represent 68%, 90%, 99% confidence level, respectively. The asterisk 
in Fig. 1 (a) stands for the best fit point. The value of Am 2 3 for each set of fig- 
ures is that of the best fit point with GIVEN value of S\. The best fit with a fixed 
value of S 1 = 0, tt/4, tt/2 is obtained for Am 2 3 = 10- 29 / 10 eV 2 , l(T 27 / 10 eV 2 , l(r 29 / 10 eV 2 , 
respectively. 

Fig. 2 (a), (b) The zenith angles dependence of (a) contained events and (b) upward going 
through muon events for Am 2 3 = 10~ 29 / 10 eV 2 , 9 2A = 50°, 9 U = 90°, 9 23 = 30° (fine 
dotted line), Am 2 3 = l(T 29 / 10 eV 2 , 9 24 = 35°, 9 M = 15°, 9 23 = 20° (dashed line; best 
fit case) Am 2 3 = l(r 27 / 10 eV 2 , 9 24 = 40°, 9 U = 0°, # 23 = 0° (coarse dotted line; best 
fit case among pure <-> v T oscillations), no oscillation case (solid line), respectively. 
81 — for all the three cases. The zenith angle dependence of the multi-GeV /i-like 
events is different for the three sets of the oscillation parameters, but that of the 
upward going muon events is almost similar for the best fit case and the pure <-> u T 
case. 

Fig. 3 (a), (b), (c) The shadowed area is the allowed region projected on the (#34, #23) 
plane for various values of Am| 3 (lCT 3 - 5 eV 2 < Am] 3 < lCT 2 eV 2 ) for each value of 
#24 = 27°,30°,---,58° and for (a) 61 = 0, (b) ^ = tt/4, (c) 5 x = tt/2, respectively, 
and the thin solid lines are boundary of the allowed region for various values of Am 2 3 . 
The solid, dashed, coarse dotted and fine dotted lines stand for the contours of c s = 
\U s i\ 2 + \U s2 \ 2 HC23C34 + S23S24S34e* <51 1 2 = 0.2,0.4,0.6,0.8, respectively. Solutions with 
c s < 0.2 exist for 40° < 6* 2 4 ^ 50° and for each value of 5± and they can have Large 
Mixing Angle solutions of the solar neutrino problem. 
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